A Single-Institution Prospective Study To Evaluate the Safety and Efficacy of Real- Time Image-Gated Spot-Scanning Proton Therapy (RGPT) for Prostate Cancer.

Purpose: In real-time image-gated spot-scanning proton therapy (RGPT), the dose distribution is distorted by gold fiducial markers placed in the prostate. Distortion can be suppressed by using small markers and more than 2 fields, but additional fields may increase the dose to organs at risk. Therefore, we conducted a prospective study to evaluate the safety and short-term clinical outcome of RGPT for prostate cancer. Methods and Materials: Based on the previously reported frequency of early adverse events (AE) and the noninferiority margin of 10%, the required number of cases was calculated to be 43 using the one-sample binomial test by the Southwest Oncology Group statistical tools with the one-sided significance level of 2.5% and the power 80%. Patients with localized prostate cancer were enrolled and 3 to 4 pure gold fiducial markers of 1.5-mm diameter were inserted in the prostate. The prescribed dose was 70 Gy(relative biologic effectiveness) in 30 fractions, and treatment was performed with 3 fields from the left, right, and the back, or 4 fields from either side of slightly anterior and posterior oblique fields. The primary endpoint was the frequency of early AE (≥grade 2) and the secondary endpoint was the biochemical relapse-free survival rate and the frequency of late AE. Results: Forty-five cases were enrolled between 2015 and 2017, and all patients completed the treatment protocol. The median follow-up period was 63.0 months. The frequency of early AE (≥grade 2) was observed in 4 cases (8.9%), therefore the noninferiority was verified. The overall 5-year biochemical relapse-free survival rate was 88.9%. As late AE, grade 2 rectal bleeding was observed in 8 cases (17.8%). Conclusions: The RGPT for prostate cancer with 1.5-mm markers and 3- or 4- fields was as safe as conventional proton therapy in early AE, and its efficacy was comparable with previous studies.

Probability of normal tissue complications for hematologic and gastrointestinal toxicity in postoperative whole pelvic radiotherapy for gynecologic malignancies using intensity-modulated proton therapy with robust optimization.

This retrospective treatment-planning study was conducted to determine whether intensity-modulated proton therapy with robust optimization (ro-IMPT) reduces the risk of acute hematologic toxicity (H-T) and acute and late gastrointestinal toxicity (GI-T) in postoperative whole pelvic radiotherapy for gynecologic malignancies when compared with three-dimensional conformal radiation therapy (3D-CRT), intensity-modulated X-ray (IMXT) and single-field optimization proton beam (SFO-PBT) therapies. All plans were created for 13 gynecologic-malignancy patients. The prescribed dose was 45 GyE in 25 fractions for 95% planning target volume in 3D-CRT, IMXT and SFO-PBT plans and for 99% clinical target volume (CTV) in ro-IMPT plans. The normal tissue complication probability (NTCP) of each toxicity was used as an in silico surrogate marker. Median estimated NTCP values for acute H-T and acute and late GI-T were 0.20, 0.94 and 0.58 × 10-1 in 3D-CRT; 0.19, 0.65 and 0.24 × 10-1 in IMXT; 0.04, 0.74 and 0.19 × 10-1 in SFO-PBT; and 0.06, 0.66 and 0.15 × 10-1 in ro-IMPT, respectively. Compared with 3D-CRT and IMXT plans, the ro-IMPT plan demonstrated significant reduction in acute H-T and late GI-T. The risk of acute GI-T in ro-IMPT plan is equivalent with IMXT plan. The ro-IMPT plan demonstrated potential clinical benefits for reducing the risk of acute H-T and late GI-T in the treatment of gynecologic malignances by reducing the dose to the bone marrow and bowel bag while maintaining adequate dose coverage to the CTV. Our results indicated that ro-IMPT may reduce acute H-T and late GI-T risk with potentially improving outcomes for postoperative gynecologic-malignancy patients with concurrent chemotherapy.

Technical note: Application of an optical hydrophone to ionoacoustic range detection in a tissue-mimicking agar phantom.

BACKGROUND: Ionoacoustics is a promising approach to reduce the range uncertainty in proton therapy. A miniature-sized optical hydrophone (OH) was used as a measuring device to detect weak ionoacoustic signals with a high signal-to-noise ratio in water. However, further development is necessary to prevent wave distortion because of nearby acoustic impedance discontinuities while detection is conducted on the patient's skin. PURPOSE: A prototype of the probe head attached to an OH was fabricated and the required dimensions were experimentally investigated using a 100-MeV proton beam from a fixed-field alternating gradient accelerator and k-Wave simulations. The beam range of the proton in a tissue-mimicking phantom was estimated by measuring γ-waves and spherical ionoacoustic waves with resonant frequency (SPIRE). METHODS: Four sizes of probe heads were fabricated from agar blocks for the OH. Using the prototype, the Î³-wave was detected at distal and lateral positions to the Bragg peak on the phantom surface for proton beams delivered at seven positions. For SPIRE, independent measurements were performed at distal on- and off-axis positions. The range positions were estimated by solving the linear equation using the sensitive matrix for the γ-wave and linear fitting of the correlation curve for SPIRE; they were compared with those measured using a film. RESULTS: The first peak of the γ-wave was undistorted with the 3 × 3 × 3-cm3 probe head used at the on-axis and 3-cm off-axis positions. The range positions estimated by the γ-wave agreed with the film-based range in the depth direction (the maximum deviation was 0.7 mm), although a 0.6-2.1 mm deviation was observed in the lateral direction. For SPIRE, the deviation was <1 mm for the two measurement positions. CONCLUSIONS: The attachment of a relatively small-sized probe head allowed the OH to measure the beam range on the phantom surface.

Study of hepatic toxicity in small liver tumors after photon or proton therapy based on factors predicting the benefits of proton.

OBJECTIVES: In a previous study of hepatic toxicity, the following three factors were identified to predict the benefits of proton beam therapy (PBT) for hepatocellular carcinomas (HCCs) with a maximum diameter of ≤5 cm and Child-pugh grade A (CP-A): number of tumors (1 vs ≥2), the location of tumors (hepatic hilum or others), and the sum of the diameters of lesions. The aim of this study is to analyze the association between these three factors and hepatic toxicity. METHODS: We retrospectively reviewed patients of CP-A treated with PBT or photon stereotactic body radiotherapy (X-ray radiotherapy, XRT) for HCC ≤5 cm. For normal liver dose, the V5, V10, V20 (volumes receiving 5, 10, and 20 Gy at least), and the mean dose was evaluated. The albumin-bilirubin (ALBI) and CP score changes from the baseline were evaluated at 3 and 6 months after treatment. RESULTS: In 89 patients (XRT: 48, PBT: 41), those with two or three (2-3) predictive factors were higher normal liver doses than with zero or one (0-1) factor. In the PBT group, the ALBI score worsened more in patients with 2-3 factors than those with 0-1 factor, at 3 months (median: 0.26 vs 0.02, p = 0.032) and at 6 months (median: 0.35 vs 0.10, p = 0.009). The ALBI score change in the XRT group and CP score change in either modality were not significantly different in the number of predictive factors. CONCLUSION: The predictive factor numbers predicted the ALBI score change in PBT but not in XRT. ADVANCES IN KNOWLEDGE: This study suggest that the number of predictive factors previously identified (0-1 vs 2-3) were significantly associated with dosimetric parameters of the normal liver in both modalities. In the proton group, the number of predictive factors was associated with a worsening ALBI score at 3 and 6 months, but these associations were not found in the photon SBRT group.

Ionoacoustic application of an optical hydrophone to detect proton beam range in water.

BACKGROUND: Proton range uncertainty has been the main factor limiting the ability of proton therapy to concentrate doses to tumors to their full potential. Ionoacoustic (IA) range verification is an approach to reducing this uncertainty by detecting thermoacoustic waves emitted from an irradiated volume immediately following a pulsed proton beam delivery; however, the signal weakness has been an obstacle to its clinical application. To increase the signal-to-noise ratio (SNR) with the conventional piezoelectric hydrophone (PH), the detector-sensitive volume needs to be large, but it could narrow the range of available beam angles and disturb real-time images obtained during beam delivery. PURPOSE: To prevent this issue, we investigated a millimeter-sized optical hydrophone (OH) that exploits the laser interferometric principle. For two types of IA waves [γ-wave emitted from the Bragg peak (BP) and a spherical IA wave with resonant frequency (SPIRE) emitted from the gold fiducial marker (GM)], comparisons were made with PH in terms of waveforms, SNR, range detection accuracy, and signal intensity robustness against the small detector misalignment, particularly for SPIRE. METHODS: A 100-MeV proton beam with a 27 ns pulse width and 4 mm beam size was produced using a fixed-field alternating gradient accelerator and was irradiated to the water phantom. The GM was set on the beam's central axis. Acrylic plates of various thicknesses, up to 12 mm, were set in front of the phantoms to shift the proton range. OH was set distal and lateral to the beam, and the range was estimated using the time-of-flight method for γ-wave and by comparing with the calibration data (SPIRE intensity versus the distance between the GM and BP) derived from an IA wave transport simulation for SPIRE. The BP dose per pulse was 0.5-0.6 Gy. To measure the variation in SPIRE amplitude against the hydrophone misalignment, the hydrophone was shifted by ± 2 mm at a maximum in lateral directions. RESULTS: Despite its small size, OH could detect γ-wave with a higher SNR than the conventional PH (diameter, 29 mm), and a single measurement was sufficient to detect the beam range with a submillimeter accuracy in water. In the SPIRE measurement, OH was far more robust against the detector misalignment than the focused PH (FPH) used in our previous study [5%/mm (OH) versus 80%/mm (FPH)], and the correlation between the measured SPIRE intensity and the distance between the GM and BP agreed well with the simulation results. However, the OH sensitivity was lower than the FPH sensitivity, and about 5.6-Gy dose was required to decrease the intensity variation among measurements to less than 10%. CONCLUSION: The miniature OH was found to detect weak IA signals produced by proton beams with a BP dose used in hypofractionated regimens. The OH sensitivity improvement at the MHz regime is worth exploring as the next step.

Technical note: Performance evaluation of volumetric imaging based on motion modeling by principal component analysis.

PURPOSE: To quantitatively evaluate the achievable performance of volumetric imaging based on lung motion modeling by principal component analysis (PCA). METHODS: In volumetric imaging based on PCA, internal deformation was represented as a linear combination of the eigenvectors derived by PCA of the deformation vector fields evaluated from patient-specific four-dimensional-computed tomography (4DCT) datasets. The volumetric image was synthesized by warping the reference CT image with a deformation vector field which was evaluated using optimal principal component coefficients (PCs). Larger PCs were hypothesized to reproduce deformations larger than those included in the original 4DCT dataset. To evaluate the reproducibility of PCA-reconstructed volumetric images synthesized to be close to the ground truth as possible, mean absolute error (MAE), structure similarity index measure (SSIM) and discrepancy of diaphragm position were evaluated using 22 4DCT datasets of nine patients. RESULTS: Mean MAE and SSIM values for the PCA-reconstructed volumetric images were approximately 80 HU and 0.88, respectively, regardless of the respiratory phase. In most test cases including the data of which motion range was exceeding that of the modeling data, the positional error of diaphragm was less than 5 mm. The results suggested that large deformations not included in the modeling 4DCT dataset could be reproduced. Furthermore, since the first PC correlated with the displacement of the diaphragm position, the first eigenvector became the dominant factor representing the respiration-associated deformations. However, other PCs did not necessarily change with the same trend as the first PC, and no correlation was observed between the coefficients. Hence, randomly allocating or sampling these PCs in expanded ranges may be applicable to reasonably generate an augmented dataset with various deformations. CONCLUSIONS: Reasonable accuracy of image synthesis comparable to those in the previous research were shown by using clinical data. These results indicate the potential of PCA-based volumetric imaging for clinical applications.

Deformed dose restoration to account for tumor deformation and position changes for adaptive proton therapy.

BACKGROUND: Online adaptation during intensity-modulated proton therapy (IMPT) can minimize the effect of inter-fractional anatomical changes, but remains challenging because of the complex workflow. One approach for fast and automated online IMPT adaptation is dose restoration, which restores the initial dose distribution on the updated anatomy. However, this method may fail in cases where tumor deformation or position changes occur. PURPOSE: To develop a fast and robust IMPT online adaptation method named "deformed dose restoration (DDR)" that can adjust for inter-fractional tumor deformation and position changes. METHODS: The DDR method comprises two steps: (1) calculation of the deformed dose distribution, and (2) restoration of the deformed dose distribution. First, the deformable image registration (DIR) between the initial clinical target volume (CTV) and the new CTV were performed to calculate the vector field. To ensure robustness for setup and range uncertainty and the ability to restore the deformed dose distribution, an expanded CTV-based registration to maintain the dose gradient outside the CTV was developed. The deformed dose distribution was obtained by applying the vector field to the initial dose distribution. Then, the voxel-by-voxel dose difference optimization was performed to calculate beam parameters that restore the deformed dose distribution on the updated anatomy. The optimization function was the sum of total dose differences and dose differences of each field to restore the initial dose overlap of each field. This method only requires target contouring, which eliminates the need for organs at risk (OARs) contouring. Six clinical cases wherein the tumor deformation and/or position changed on repeated CTs were selected. DDR feasibility was evaluated by comparing the results with those from three other strategies, namely, not adapted (continuing the initial plan), adapted by previous dose restoration, and fully optimized. RESULTS: In all cases, continuing the initial plan was largely distorted on the repeated CTs and the dose-volume histogram (DVH) metrics for the target were reduced due to the tumor deformation or position changes. On the other hand, DDR improved DVH metrics for the target to the same level as the initial dose distribution. Dose increase was seen for some OARs because tumor growth had reduced the relative distance between CTVs and OARs. Robustness evaluation for setup and range uncertainty (3 mm/3.5%) showed that deviation in DVH-bandwidth for CTV D95% from the initial plan was 0.4% ± 0.5% (Mean ± S.D.) for DDR. The calculation time was 8.1 ± 6.4 min. CONCLUSIONS: An online adaptation algorithm was developed that improved the treatment quality for inter-fractional anatomical changes and retained robustness for intra-fractional setup and range uncertainty. The main advantage of this method is that it only requires target contouring alone and saves the time for OARs contouring. The fast and robust adaptation method for tumor deformation and position changes described here can reduce the need for offline adaptation and improve treatment efficiency.

Hyperfractionated intensity-modulated proton therapy for pharyngeal cancer with variable relative biological effectiveness: A simulation study.

BACKGROUND: The relative biological effectiveness (RBE) of proton is considered to be dependent on biological parameters and fractional dose. While hyperfractionated photon therapy was effective in the treatment of patients with head and neck cancers, its effect in intensity-modulated proton therapy (IMPT) under the variable RBE has not been investigated in detail. PURPOSE: To study the effect of variable RBE on hyperfractionated IMPT for the treatment of pharyngeal cancer. We investigated the biologically effective dose (BED) to determine the theoretical effective hyperfractionated schedule. METHODS: The treatment plans of three pharyngeal cancer patients were used to define the ΔBED for the clinical target volume (CTV) and soft tissue (acute and late reaction) as the difference between the BED for the altered schedule with variable RBE and conventional schedule with constant RBE. The ΔBED with several combinations of parameters (treatment days, number of fractions, and prescribed dose) was comprehensively calculated. Of the candidate schedules, the one that commonly gave a higher ΔBED for CTV was selected as the resultant schedule. The BED volume histogram was used to compare the influence of variable RBE and fractionation. RESULTS: In the conventional schedule, compared with the constant RBE, the variable RBE resulted in a mean 2.6 and 2.7 Gy reduction of BEDmean for the CTV and soft tissue (acute reaction) of the three plans, respectively. Moreover, the BEDmean for soft tissue (late reaction) increased by 7.4 Gy, indicating a potential risk of increased RBE. Comprehensive calculation of the ΔBED resulted in the hyperfractionated schedule of 80.52 Gy (RBE = 1.1)/66 fractions in 6.5 weeks. When variable RBE was used, compared with the conventional schedule, the hyperfractionated schedule increased the BEDmean for CTV by 7.6 Gy; however, this was associated with a 7.8 Gy increase for soft tissue (acute reaction). The BEDmean for soft tissue (late reaction) decreased by 2.4 Gy. CONCLUSION: The results indicated a potential effect of the variable RBE on IMPT for pharyngeal cancer but with the possibility that hyperfractionation could outweigh this effect. Although biological uncertainties require conservative use of the resultant schedule, hyperfractionation is expected to be an effective strategy in IMPT for pharyngeal cancer.

A study on predicting cases that would benefit from proton beam therapy in primary liver tumors of less than or equal to 5 cm based on the estimated incidence of hepatic toxicity.

Background: For small primary liver tumors, favorable outcomes have been reported with both of proton beam therapy (PBT) and X-ray therapy (XRT). However, no clear criteria have been proposed in the cases for which and when of PBT or XRT has to be used. The aim of this study is to investigate cases that would benefit from PBT based on the predicted rate of hepatic toxicity. Materials and methods: Eligible patients were those who underwent PBT for primary liver tumors with a maximum diameter of ≤ 5 cm and Child-Pugh grade A (n = 40). To compare the PBT-plan, the treatment plan using volumetric modulated arc therapy was generated as the XRT-plan. The rate of predicted hepatic toxicity was estimated using five normal tissue complication probability (NTCP) models with three different endpoints. The differences in NTCP values (ΔNTCP) were calculated to determine the relative advantage of PBT. Factors predicting benefits of PBT were analyzed by logistic regression analysis. Results: From the dose-volume histogram comparisons, an advantage of PBT was found in sparing of the normal liver receiving low doses. The factors predicting the benefit of PBT differed depending on the selected NTCP model. From the five models, the total tumor diameter (sum of the target tumors), location (hepatic hilum vs other), and number of tumors (1 vs 2) were significant factors. Conclusions: From the radiation-related hepatic toxicity, factors were identified to predict benefits of PBT in primary liver tumors with Child-Pugh grade A, with the maximum tumor diameter of ≤ 5 cm.

A Consistent Protocol Reveals a Large Heterogeneity in the Biological Effectiveness of Proton and Carbon-Ion Beams for Various Sarcoma and Normal-Tissue-Derived Cell Lines.

This study investigated variations in the relative biological effectiveness (RBE) values among various sarcoma and normal-tissue-derived cell lines (normal cell line) in proton beam and carbon-ion irradiations. We used a consistent protocol that specified the timing of irradiation after plating cells and detailed the colony formation assay. We examined the cell type dependence of RBE for proton beam and carbon-ion irradiations using four human sarcoma cell lines (MG63 osteosarcoma, HT1080 fibrosarcoma, SW872 liposarcoma, and SW1353 chondrosarcoma) and three normal cell lines (HDF human dermal fibroblast, hTERT-HME1 mammary gland, and NuLi-1 bronchus epithelium). The cells were irradiated with gamma rays, proton beams at the center of the spread-out Bragg peak, or carbon-ion beams at 54.4 keV/µm linear energy transfer. In all sarcoma and normal cell lines, the average RBE values in proton beam and carbon-ion irradiations were 1.08 ± 0.11 and 2.08 ± 0.36, which were consistent with the values of 1.1 and 2.13 used in current treatment planning systems, respectively. Up to 34% difference in the RBE of the proton beam was observed between MG63 and HT1080. Similarly, a 32% difference in the RBE of the carbon-ion beam was observed between SW872 and the other sarcoma cell lines. In proton beam irradiation, normal cell lines had less variation in RBE values (within 10%), whereas in carbon-ion irradiation, RBE values differed by up to 48% between hTERT-HME1 and NuLi-1. Our results suggest that specific dose evaluations for tumor and normal tissues are necessary for treatment planning in both proton and carbon-ion therapies.

An initial systematic study of the linear energy transfer distributions of a proton beam under a transverse magnetic field.

PURPOSE: To evaluate the biological effectiveness of magnetic resonance (MR)-guided proton beam therapy, comprehensively characterizing the dose and dose-averaged linear energy transfer (LETd ) distributions under a magnetic field is necessary. Although detailed analysis has characterized curved beam paths and distorted dose distributions, the impact of a magnetic field on LETd should also be explored to determine the proton relative biological effectiveness (RBE). Hence, this initial study aims to present a basic analysis of LETd distributions in the presence of a magnetic field using Monte Carlo simulation (MCS). METHODS: Geant4 MCS (version 10.1.p01) was performed to calculate the LETd distribution of proton beams. The incident beam energies were set to 70.2, 140.8, and 220 MeV, and both zero- and finite-emittance pencil beams as well as scanned field were simulated. A transverse magnetic field of 0-3 T was applied within a water phantom placed at the isocenter, and the three-dimensional dose and LETd distributions in the phantom were calculated. Then, the depth profiles of LETd along the curved trajectory and the lateral LETd profile at the Bragg peak (BP) depth were analyzed under changing energies and magnetic fields. In addition, for zero- and finite-emittance beams, the correlation of the lateral asymmetries between the dose and LETd distributions were analyzed. Finally, spread-out Bragg peak (SOBP) fields were simulated to assess the depth-dependent asymmetry of the LETd distributions. RESULTS: A transverse magnetic field distorted the lateral LETd distribution of a pencil beam at close to the BP, and the magnitude of the distortion at the BP increased for higher energy beams and larger magnetic fields. For a zero-emittance beam, the differences in LETd between the left and right D20 positions were relatively large; the difference in LETd was 1.5 and 2.3 keV/µm at 140.8 and 220 MeV, respectively, at a magnetic field of 1.5 T. These asymmetries were pronounced at positions where the dose asymmetries were large. The size of the asymmetry was less substantial for a finite-emittance beam and even less for a scanned field. However, a 1.5-keV/µm difference still remained between the left and right D20 positions of a scanned field penumbra for a 220 MeV beam under the same magnetic field. For the SOBP field, it was found that the distal region of SOBP had the highest LETd distortions, followed by the central and proximal regions for the middle-sized SOBP (5 × 5 × 5 cm3 ), whereas the degree of LETd distortion did not vary much with depth for the 10 × 10 × 10-cm3 SOBP field. CONCLUSION: Our results indicate that not only the dose but also LETd distortions should be considered to accurately evaluate the biological effectiveness of MR-guided proton beam therapy.

Particle therapy using protons or carbon ions for cancer patients with cardiac implantable electronic devices (CIED): a retrospective multi-institutional study.

PURPOSE: To evaluate the outcomes of particle therapy in cancer patients with cardiac implantable electronic devices (CIEDs). MATERIALS AND METHODS: From April 2001 to March 2013, 19,585 patients were treated with proton beam therapy (PBT) or carbon ion therapy (CIT) at 8 institutions. Of these, 69 patients (0.4%, PBT 46, CIT 22, and PBT + CIT 1) with CIEDs (64 pacemakers, 4 implantable cardioverter defibrillators, and 1 with a cardiac resynchronization therapy defibrillator) were retrospectively reviewed. All the patients with CIEDs in this study were treated with the passive scattering type of particle beam therapy. RESULTS: Six (13%) of the 47 PBT patients, and none of the 23 CIT patients experienced CIED malfunctions (p = 0.105). Electrical resets (7) and over-sensing (3) occurred transiently in 6 patients. The distance between the edge of the irradiation field and the CIED was not associated with the incidence of malfunctions in 20 patients with lung cancer. A larger field size had a higher event rate but the test to evaluate trends as not statistically significant (p = 0.196). CONCLUSION: Differences in the frequency of occurrence of device malfunctions for patients treated with PBT and patients treated with CIT did not reach statistical significance. The present study can be regarded as a benchmark study about the incidence of malfunctioning of CIED in passive scattering particle beam therapy and can be used as a reference for active scanning particle beam therapy.

Impact of a spatially dependent dose delivery time structure on the biological effectiveness of scanning proton therapy.

PURPOSE: In the scanning beam delivery of protons, different portions of the target are irradiated with different linear energy transfer protons with various time intervals and irradiation times. This research aimed to evaluate the spatially dependent biological effectiveness of protracted irradiation in scanning proton therapy. METHODS: One and two parallel opposed fields plans were created in water phantom with the prescribed dose of 2 Gy. Three scenarios (instantaneous, continuous, and layered scans) were used with the corresponding beam delivery models. The biological dose (physical dose × relative biological effectiveness) was calculated using the linear quadratic model and the theory of dual radiation action to quantitatively evaluate the dose delivery time effect. In addition, simulations using clinical plans (postoperative seminoma and prostate tumor cases) were conducted to assess the impact of the effects on the dose volume histogram parameters and homogeneity coefficient (HC) in targets. RESULTS: In a single-field plan of water phantom, when the treatment time was 19 min, the layered-scan scenario showed a decrease of <0.2% (almost 3.3%) in the biological dose from the plan on the distal (proximal) side because of the high (low) dose rate. This is in contrast to the continuous scenario, where the biological dose was almost uniformly decreased over the target by approximately 3.3%. The simulation with clinical geometry showed that the decrease rates in D99% were 0.9% and 1.5% for every 10 min of treatment time prolongation for postoperative seminoma and prostate tumor cases, respectively, whereas the increase rates in HC were 0.7% and 0.2%. CONCLUSIONS: In protracted irradiation in scanning proton therapy, the spatially dependent dose delivery time structure in scanning beam delivery can be an important factor for accurate evaluation of biological effectiveness.

A treatment planning study of urethra-sparing intensity-modulated proton therapy for localized prostate cancer.

BACKGROUND AND PURPOSE: Urethra-sparing radiation therapy for localized prostate cancer can reduce the risk of radiation-induced genitourinary toxicity by intentionally underdosing the periurethral transitional zone. We aimed to compare the clinical impact of a urethra-sparing intensity-modulated proton therapy (US-IMPT) plan with that of conventional clinical plans without urethral dose reduction. MATERIALS AND METHODS: This study included 13 patients who had undergone proton beam therapy. The prescribed dose was 63 GyE in 21 fractions for 99% of the clinical target volume. To compare the clinical impact of the US-IMPT plan with that of the conventional clinical plan, tumor control probability (TCP) and normal tissue complication probability (NTCP) were calculated with a generalized equivalent uniform dose-based Lyman-Kutcher model using dose volume histograms. The endpoints of these model parameters for the rectum, bladder, and urethra were fistula, contraction, and urethral stricture, respectively. RESULTS: The mean NTCP value for the urethra in US-IMPT was significantly lower than that in the conventional clinical plan (0.6% vs. 1.2%, p < 0.05). There were no statistically significant differences between the conventional and US-IMPT plans regarding the mean minimum dose for the urethra with a 3-mm margin, TCP value, and NTCP value for the rectum and bladder. Additionally, the target dose coverage of all plans in the robustness analysis was within the clinically acceptable range. CONCLUSIONS: Compared with the conventional clinically applied plans, US-IMPT plans have potential clinical advantages and may reduce the risk of genitourinary toxicities, while maintaining the same TCP and NTCP in the rectum and bladder.

First experimental results of gated proton imaging using x-ray fluoroscopy to detect a fiducial marker.

Increasing numbers of proton imaging research studies are being conducted for accurate proton range determination in proton therapy treatment planning. However, there is no proton imaging system that deals with motion artifacts. In this study, a gated proton imaging system was developed and the first experimental results of proton radiography (pRG) were obtained for a moving object without motion artifacts. A motion management system using dual x-ray fluoroscopy for detecting a spherical gold fiducial marker was introduced and the proton beam was gated in accordance with the motion of the object. To demonstrate the performance of the gated proton imaging system, gated pRG images of a moving phantom were acquired experimentally, and the motion artifacts clearly were diminished. Also, the factors causing image deteriorations were evaluated focusing on the new gating system developed here, and the main factor was identified as the latency (with a maximum value of 93 ms) between the ideal gating signal according to the actual marker position and the actual gating signal. The possible deterioration due to the latency of the proton imaging system and proton beam irradiation was small owing to appropriate setting of the time structure.

Real-time CT image generation based on voxel-by-voxel modeling of internal deformation by utilizing the displacement of fiducial markers.

PURPOSE: To show the feasibility of real-time CT image generation technique utilizing internal fiducial markers that facilitate the evaluation of internal deformation. METHODS: In the proposed method, a linear regression model that can derive internal deformation from the displacement of fiducial markers is built for each voxel in the training process before the treatment session. Marker displacement and internal deformation are derived from the four-dimensional computed tomography (4DCT) dataset. In the treatment session, the three-dimensional deformation vector field is derived according to the marker displacement, which is monitored by the real-time imaging system. The whole CT image can be synthesized by deforming the reference CT image with a deformation vector field in real-time. To show the feasibility of the technique, image synthesis accuracy and tumor localization accuracy were evaluated using the dataset generated by extended NURBS-Based Cardiac-Torso (XCAT) phantom and clinical 4DCT datasets from six patients, containing 10 CT datasets each. In the validation with XCAT phantom, motion range of the tumor in training data and validation data were about 10 and 15 mm, respectively, so as to simulate motion variation between 4DCT acquisition and treatment session. In the validation with patient 4DCT dataset, eight CT datasets from the 4DCT dataset were used in the training process. Two excluded inhale CT datasets can be regarded as the datasets with large deformations more than training dataset. CT images were generated for each respiratory phase using the corresponding marker displacement. Root mean squared error (RMSE), normalized RMSE (NRMSE), and structural similarity index measure (SSIM) between the original CT images and the synthesized CT images were evaluated as the quantitative indices of the accuracy of image synthesis. The accuracy of tumor localization was also evaluated. RESULTS: In the validation with XCAT phantom, the mean NRMSE, SSIM, and three-dimensional tumor localization error were 7.5 ± 1.1%, 0.95 ± 0.02, and 0.4 ± 0.3 mm, respectively. In the validation with patient 4DCT dataset, the mean RMSE, NRMSE, SSIM, and three-dimensional tumor localization error in six patients were 73.7 ± 19.6 HU, 9.2 ± 2.6%, 0.88 ± 0.04, and 0.8 ± 0.6 mm, respectively. These results suggest that the accuracy of the proposed technique is adequate when the respiratory motion is within the range of the training dataset. In the evaluation with a marker displacement larger than that of the training dataset, the mean RMSE, NRMSE, and tumor localization error were about 100 HU, 13%, and <2.0 mm, respectively, except for one case having large motion variation. The performance of the proposed method was similar to those of previous studies. Processing time to generate the volumetric image was <100 ms. CONCLUSION: We have shown the feasibility of the real-time CT image generation technique for volumetric imaging.

Prediction of target position from multiple fiducial markers by partial least squares regression in real-time tumor-tracking radiation therapy.

The purpose of this work is to show the usefulness of a prediction method of tumor location based on partial least squares regression (PLSR) using multiple fiducial markers. The trajectory data of respiratory motion of four internal fiducial markers inserted in lungs were used for the analysis. The position of one of the four markers was assumed to be the tumor position and was predicted by other three fiducial markers. Regression coefficients for prediction of the position of the tumor-assumed marker from the fiducial markers' positions is derived by PLSR. The tracking error and the gating error were evaluated assuming two possible variations. First, the variation of the position definition of the tumor and the markers on treatment planning computed tomograhy (CT) images. Second, the intra-fractional anatomical variation which leads the distance change between the tumor and markers during the course of treatment. For comparison, rigid predictions and ordinally multiple linear regression (MLR) predictions were also evaluated. The tracking and gating errors of PLSR prediction were smaller than those of other prediction methods. Ninety-fifth percentile of tracking/gating error in all trials were 3.7/4.1 mm, respectively in PLSR prediction for superior-inferior direction. The results suggested that PLSR prediction was robust to variations, and clinically applicable accuracy could be achievable for targeting tumors.

Technical Note: Range verification of pulsed proton beams from fixed-field alternating gradient accelerator by means of time-of-flight measurement of ionoacoustic waves.

PURPOSE: Ionoacoustics is one of the promising approaches to verify the beam range in proton therapy. However, the weakness of the wave signal remains a main hindrance to its application in clinics. Here we studied the potential use of a fixed-field alternating gradient accelerator (FFA), one of the accelerator candidates for future proton therapy. For such end, magnitude of the pressure wave and range accuracy achieved by the short-pulsed beam of FFA were assessed, using both simulation and experimental procedure. METHODS: A 100 MeV proton beam from the FFA was applied on a water phantom, through the acrylic wall. The beam range measured by the Bragg peak (BP)-ionization chamber (BPC) was 77.6 mm, while the maximum dose at BP was estimated to be 0.35 Gy/pulse. A hydrophone was placed 20 mm downstream of the BP, and signals were amplified and stored by a digital oscilloscope, averaged, and low-pass filtered. Time-of-flight (TOF) and two relative TOF values were analyzed in order to determine the beam range. Furthermore, an acoustic wave transport simulation was conducted to estimate the amplitude of the pressure waves. RESULTS: The range calculated when using two relative TOF was 78.16 ± 0.01 and 78.14 ± 0.01 mm, respectively, both values being coherent with the range measured by the BPC (the difference was 0.5-0.6 mm). In contrast, utilizing the direct TOF resulted in a range error of 1.8 mm. Fivefold and 50-fold averaging were required to suppress the range variation to below 1 mm for TOF and relative TOF measures, respectively. The simulation suggested the magnitude of pressure wave at the detector exceeded 7 Pascal. CONCLUSION: A submillimeter range accuracy was attained with a pulsed beam of about 21 ns from an FFA, at a clinical energy using relative TOF. To precisely quantify the range with a single TOF measurement, subsequent improvement in the measuring system is required.

Validation of dose distribution for liver tumors treated with real-time-image gated spot-scanning proton therapy by log data based dose reconstruction.

In spot scanning proton therapy (SSPT), the spot position relative to the target may fluctuate through tumor motion even when gating the radiation by utilizing a fiducial marker. We have established a procedure that evaluates the delivered dose distribution by utilizing log data on tumor motion and spot information. The purpose of this study is to show the reliability of the dose distributions for liver tumors treated with real-time-image gated SSPT (RGPT). In the evaluation procedure, the delivered spot information and the marker position are synchronized on the basis of log data on the timing of the spot irradiation and fluoroscopic X-ray irradiation. Then a treatment planning system reconstructs the delivered dose distribution. Dose distributions accumulated for all fractions were reconstructed for eight liver cases. The log data were acquired in all 168 fractions for all eight cases. The evaluation was performed for the values of maximum dose, minimum dose, D99, and D5-D95 for the clinical target volumes (CTVs) and mean liver dose (MLD) scaled by the prescribed dose. These dosimetric parameters were statistically compared between the planned dose distribution and the reconstructed dose distribution. The mean difference of the maximum dose was 1.3% (95% confidence interval [CI]: 0.6%-2.1%). Regarding the minimum dose, the mean difference was 0.1% (95% CI: -0.5%-0.7%). The mean differences of D99, D5-D95 and MLD were below 1%. The reliability of dose distributions for liver tumors treated with RGPT-SSPT was shown by the evaluation of the accumulated dose distributions.

Potential benefits of adaptive intensity-modulated proton therapy in nasopharyngeal carcinomas.

PURPOSE: To investigate potential advantages of adaptive intensity-modulated proton beam therapy (A-IMPT) by comparing it to adaptive intensity-modulated X-ray therapy (A-IMXT) for nasopharyngeal carcinomas (NPC). METHODS: Ten patients with NPC treated with A-IMXT (step and shoot approach) and concomitant chemotherapy between 2014 and 2016 were selected. In the actual treatment, 46 Gy in 23 fractions (46Gy/23Fx.) was prescribed using the initial plan and 24Gy/12Fx was prescribed using an adapted plan thereafter. New treatment planning of A-IMPT was made for the same patients using equivalent dose fractionation schedule and dose constraints. The dose volume statistics based on deformable images and dose accumulation was used in the comparison of A-IMXT with A-IMPT. RESULTS: The means of the Dmean of the right parotid gland (P < 0.001), right TM joint (P < 0.001), left TM joint (P < 0.001), oral cavity (P < 0.001), supraglottic larynx (P = 0.001), glottic larynx (P < 0.001), , middle PCM (P = 0.0371), interior PCM (P < 0.001), cricopharyngeal muscle (P = 0.03643), and thyroid gland (P = 0.00216), in A-IMPT are lower than those of A-IMXT, with statistical significance. The means of, D0.03cc , and Dmean of each sub portion of auditory apparatus and D30% for Eustachian tube and D0.5cc for mastoid volume in A-IMPT are significantly lower than those of A-IMXT. The mean doses to the oral cavity, supraglottic larynx, and glottic larynx were all reduced by more than 20 Gy (RBE = 1.1). CONCLUSIONS: An adaptive approach is suggested to enhance the potential benefit of IMPT compared to IMXT to reduce adverse effects for patients with NPC.